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Big Bang nucleosynthesis (BBN) theory predicts the primordial abundances of the light elements ? H (referred 
to as deuterium, or D for short), 3 He, 4H e, and "Li produced in the early universe. Among these, deuterium, the 
first nuclide produced by BBN, is a key primordial material for subsequent reactions. To date, the uncertainty 
in predicted deuterium abundance (D/H) remains larger than the observational precision. In this study, the 
Monte Carlo simulation code PRIMAT was used to investigate the sensitivity of 11 important BBN reactions to 
deuterium abundance. We found that the reaction rate uncertainties of the four reactions d(d, n)? He, d(d, p)t, 
d(p, y)? He, and p(n, y)d had the largest influence on the calculated D/H uncertainty. Currently, the calculated 
D/H uncertainty cannot reach observational precision even with the recent LUNA precise d(p, y)’ H e rate. 
From the nuclear physics aspect, there is still room to largely reduce the reaction-rate uncertainties; hence, 
further measurements of the important reactions involved in BBN are still necessary. A photodisintegration 
experiment will be conducted at the Shanghai Laser Electron Gamma Source (SLEGS) Facility to precisely 


study the deuterium production reaction of p(n, y)d. 


Keywords: Big Bang Nucleosynthesis; Abundance of deuterium; Reaction cross section; Reaction rate; Monte Carlo method. 


I. INTRODUCTION 


The hot Big Bang theory, first proposed in 1946 by 
Gamow [1], is now the most widely accepted cosmological 
model of the universe, where it has expanded from a very 
high- density state dominated by radiation. This theory has 
been confirmed by the observations of the cosmic microwave 
background [2-4], the expansion of the universe, and good 
global agreement between the predictions and observations 
of the primordial abundances of the lightest elements in na- 
ture: hydrogen, helium, and lithium. According to the Big 
Bang theory, the universe began with a fireball approximately 
13.8 billion years ago. Following inflation and cooling, pri- 
mordial big-bang nucleosynthesis(BBN) began when the uni- 
verse was approximately 3 min old (when the temperature 
was reduced to approximately 1 GK, i.e., particle energy 
E x kT = 0.1 MeV) and ended less than half an hour later 
when the nuclear reactions were quenched by the low temper- 
ature and density conditions in the expanding universe. Only 
the lightest nuclides were synthesized in appreciable quanti- 
ties through BBN: approximately 75% 'H and 25% “He, with 
small amounts of 2H, “He, and “Li. These relics provide us 
with a unique window into the early universe. More compre- 
hensive reviews on BBN can be found in the literature. [5-7]. 

In general, the primordial abundances of 7H (referred to as 
D) and He inferred from observational data agree with pre- 
dictions, except for the lithium problem [8, 9]. Deuterium, a 
fragile isotope, was destroyed after the BBN. Its most primi- 
tive abundance was determined from observing cosmological 
clouds at a high redshift on the line of sight of distant quasars. 
Recently, the precision of deuterium observations in cosmo- 
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logical clouds has dramatically improved, reaching an accu- 
racy of 1.19% for primordial deuterium abundance (D/H), 
that is, D/H = (2.527+0.030)x 1075 [10]. Therefore, the nu- 
clear cross-section data relevant to the deuterium involved in 
the BBN network needs to be known with similar precision 
to constrain further the cosmological parameter, that is, the 
cosmic baryon density. In this study, we focused on BBN 
deuterium abundance. 

In BBN, deuterium synthesis is the first to occur, and the 
accumulation of primitive deuterium significantly affects the 
rate of subsequent reactions. Therefore, accurate determina- 
tion of the abundance of deuterium in BBN calculations is im- 
portant. Eextensive studies relevant to BBN deuterium abun- 
dance have been conducted [11-13] but only a few mention 
the impact of the uncertainties of the relevant reaction rates 
on D/H uncertainty. In this study, we used the Monte Carlo 
simulation code PRIMAT [14] to study the D/H uncertainty 
within BBN reaction networks and demonstrated the reaction 
rate sensitivity for some major reactions. In addition, the un- 
certainties of both the nuclear physics input and cosmological 
parameters were proposed to satisfy the accuracy of the ob- 
served deuterium abundance. 


Il. BBN MODEL & MONTE CARLO METHOD 


For a standard BBN model, the evolution of the nucleosyn- 
thetic abundance can be obtained by solving the following 
system of differential equations [15]: 


dY; 
dt dN 
j,k,l 


yey 
F$ G J 
Vij Ni! N;! 
(1) 


where N; represents the mass number of the corresponding 
nuclide, Y; is the abundance of the corresponding nuclide, and 
T denotes the reaction rate, which is usually obtained using 
the following formula [9, 16, 17]: 
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At the same time , the standard Big Bang model is also a 
single-parameter model, with only the parameter 7) used, that 
is, the ratio of the number of baryons to the number of pho- 
tons; different 7 values correspond to different nuclide abun- 
dance evolution curves or frozen abundances. In the study of 
BBN, we usually obtain the corresponding abundance value 
by entering the observed 7 into the model as a parameter and 
then comparing it with the observed abundance. To perform 
such calculations, input parameters such as nuclear reaction 
rates, element abundances, and baryon density are required. 
However, because reaction network models incorporate an in- 
creasing number of nuclides, conventional numerical calcu- 
lations have become increasingly computationally expensive 
when computing uncertainties. Furthermore, according to the 
research conducted by Longland et al. [18], the traditional 
method of error propagation fails to incorporate the statistical 
significance of errors. In this study, we adopted the Monte 
Carlo sampling technique to calculate the abundance of the 
targeted nuclide. The abundance uncertainty can be directly 
determined by generating a data distribution through the sam- 
pling process. 

The Monte Carlo method is a statistical method based on 
the large number theorem. Converting the solutions of math- 
ematical problems into random samples can significantly re- 
duce the difficulty in solving complex models. Because of 
the sampling of specific physical quantities following a cer- 
tain distribution (such as Poisson, Gaussian, and log-normal 
distributions), it also preserves the statistical significance of 
the physical quantities and their associated uncertainties and 
applies them realistically in the calculation results. 

In this study, we used the PRIMAT code [14] for BBN 
calculations. PRIMAT is mainly divided into the following 
parts: first, it needs to determine some cosmological param- 
eters, such as the parameter a(t) of the Friedmanan-Lemaitre 
(FL) spacetime, which in the code is obtained from the ther- 
modynamics of the plasma, and its change with temperature 
is obtained by numerical inversion; second, it calculates the 
effects of some weak interactions and stores the relationship 
between them and the temperature in a hard disk; finally, a 
reaction network model is established to calculate the nuclear 
reaction processes of >10 GK, 10—1.25 GK, and <1.25 GK 
according to the temperature change. During the entire pro- 
cess, the code uses random sampling of the distribution of rel- 
evant parameters to obtain the distribution function of abun- 
dance and then determines the uncertainty of the abundances. 

In particular, the main focus of BBN is nuclear physics in- 
put quantities, specifically the reaction rates of the relevant 
reactions. The PRIMAT code assumes that the reaction rate 
follows a lognormal distribution [14]: 


1 1 _nz-p)? 
f(x) = — -e z (3) 
ov 2r T 
where x denotes the reaction rate, that is, x = Na(ov), u 


and o are parameters of the lognormal distributions. For this 
distribution, the corresponding low, median, and high rates 
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are expressed as follows: 


= -—o = a +o 
Ziow = el ’ Tmed = eM, Thigh = e” (4) 
where Zjow, Zmed and high represent the reaction rates with 


probabilities of 16 


I. MODEL PARAMETERS & SIMULATION 


The main parameters in the PRIMAT code that can be ad- 


122 justed are the reaction rates and the associated uncertainties 
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of the 11 reactions of primary importance [19] involved in 
the BBN network to evaluate their impact on deuterium abun- 
dance. The reaction rates were adopted from the default rates 
used in the PRIMAT code, and the data sources are listed in 
Table 1. In addition, the default values set in the PRIMAT 
code were adopted, that is, the cosmological parameter was 
Qh? = (0.02225+0.00016) [4] for the baryon-to-photon ra- 
tio, and the neutron lifetime was Tn = (879.5-£0.8) s [20]. No- 
tably, in the original PRIMAT calculations [14], a calculated 
D/H uncertainty of 1.46% was obtained using the default re- 
action rates and their uncertainties, but without considering 
the uncertainties of the cosmological parameter 9h? at that 
time. In our calculations, a similar D/H uncertainty of 1.47% 
was obtained under the same conditions, confirming the cor- 
rectness of our calculations. 


Table 1. Data source of the default reaction rates used in the PRI- 
MAT code. Other relevant results are not included in this calcula- 
tion * 


Reaction Reference 
p(n, y)d Ando et al. [21] 
d(p, y)’ He Iliadis et al. [22] 
d(d,n)* He Gómez Inesta et al. [23] 
d(d, p)t Gómez Inesta et al. [23] 
3He(n, p)t Descouvemont et al. [24] 
t(d,n)‘He Descouvemont et al. [24] 
3 He(d, p)*He Descouvemont et al. [24] 
3 He(a, 7)" Be Tliadis et al. [22] 
tla, y) Li Descouvemont et al. [24] 
T Be(n, p)’ Li Descouvemont et al. [24] 


"Li(p, a)“ He Descouvemont et al. [24] 


a Viviani et al. [25] theoretically studied the d(d,p) and d(d,n) reactions, 
and they only plotted the astrophysical S factors and uncertainties, while 
their uncertainties were just a rough estimate. No new reaction rate and 
uncertainty were provided in that paper; Tumino et al. [26] and Pizzone et 
al. [27] reported the Trojan-Horse-Method (THM) results on the d(d,p) 
and d(d,n) reactions and gave a lower and upper limits of 4.5% and 5.0%, 
respectively. However, their data were not included in Ref. [23], which 
was adopted in the PRIMAT code. This (4.5-5.0)% error is much larger 
than the adopted one, and thus can be covered by the present sensitivity 
study. Furthermore, the full Bayesian analysis for these two reactions 
including the new THM data might also make a shift in the centroid 
values, which is beyond the scope of this study. 


First, we studied the effect of uncertainties in the 11 reac- 
tion rates, Q,h? and Tn on the calculated D/H uncertainties. 
The results are listed in Table 2. This shows that the uncer- 
tainty in the neutron lifetime Ar, negligibly impacts the D/H 
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Table 2. Sensitivity of uncertainties in Qyh?, Tn and 11 reaction 
rates on the D/H uncertainty. The results for the “He abundance 
Yp (i.e., mass fraction of 4“ He) are also listed. Here “default” de- 
notes using the default value in the PRIMAT code [14], i.e., Ah? 
= 0.00016, At, = 0.8 s, and the default reaction-rate uncertainties 
(ARates); “w/o” denotes without considering the corresponding un- 
certainty (i.e., set to zero). 


AQgh Atn ARates D/H uncer. (%) Yp uncer. (%) 
w/o w/o w/o 2.460(0.05%) 0.247(0.001%) 
default w/o w/o 2.458(1.1998%) 0.247(0.03%) 
w/o default w/o 2.460(0.06%) 0.247(0.07%) 
w/o w/o default 2.460(1.50%) 0.247(0.01%) 
w/o default default 2.460(1.47%) 0.247(0.07%) 
default w/o default 2.459(1.87%) 0.247(0.03%) 
default default w/o 2.458(1.1995%)  0.247(0.08%) 
default default default 2.459(1.83%) 0.247(0.07%) 


uncertainty, whereas the uncertainties in the reaction rates and 
Qh? have a much larger impact. Conversely, At, had a sig- 
nificant impact on the He abundance. The impact of the neu- 
tron lifetime on Big Bang nucleosynthesis was studied in Ref. 
[28] and is not discussed here. 

In the following, we investigate the sensitivity of the 
reaction-rate uncertainties to the primordial deuterium abun- 
dance (D/H) uncertainties by considering the uncertainties of 
Qah? and Tn simultaneously. The default reaction rate un- 
certainties adopted in PRIMAT for the four most important 
reactions are presented in Table 3. 


Table 3. Rates uncertainties for the four most important reactions of 
deuterium abundance. 


T(GK) p(n, 7)[21] d(p,7)[22] d(d, n)[23] _d(d, p)[23] 
0.001 0.45% 3.70% 1.10% 1.10% 
0.005 0.45% 3.70% 1.10% 1.10% 
0.010 0.45% 3.70% 1.10% 1.10% 
0.050 0.46% 3.70% 1.10% 1.10% 
0.100 0.46% 3.70% 1.10% 1.10% 
0.500 0.44% 3.70% 1.10% 1.10% 
1.000 0.44% 3.70% 1.10% 1.10% 
1.500 0.46% 3.70% 1.10% 1.10% 
2.000 0.49% 3.70% 1.10% 1.10% 
3.000 0.56% 3.70% 1.40% 1.40% 
4.000 0.61% 3.70% 1.50% 1.50% 
5.000 0.64% 4.00% 1.60% 1.60% 
6.000 0.67% 4.20% 1.70% 1.70% 
7.000 0.69% 4.40% 1.80% 1.80% 
8.000 0.71% 4.60% 1.80% 1.80% 
9.000 0.72% 4.70% 1.80% 1.80% 
10.000 0.74% 4.90% 1.80% 1.90% 


To examine the sensitivity, we multiplied the default uncer- 
tainty of just one reaction rate by arbitrary factors of M.F. = 
0.1, 0.8, 0.9, 1, 2, 4 and 6, respectively, and set the remaining 
10 reactions with their default rate uncertainties. We found 
that the d(p, y)?” He, d(d, n)? He, d(d, p)t and p(n, y)d reac- 
tions are the four major reactions that have the largest influ- 
ence on the D/H uncertainty, compared to the remaining reac- 
tions. The calculated results are listed in Table 4 and shown 


in Fig. 1. 

For multiplying factors of less than one (i.e., reducing the 
current rate uncertainty), the decreasing trend for the four re- 
actions of interest is shown in the inset of Fig. 1. This demon- 
strates that the calculated D/H uncertainties can reach 1.41%, 
1.76%, 1.77%, and 1.82% when the reaction rate uncertain- 
ties of d(p, y)? He, d(d, n)? He, d(d, p)t and p(n, y)d are re- 
duced by a factor of 10 (i.e., M.F. = 0.1, nearly no rate uncer- 
tainty). In general, the current default d(p, y)’ He rate uncer- 
tainty (i.e., 3.7%) dominates the calculated D/H uncertainty 
because the default rate uncertainties adopted for the latter 
three reactions are already quite small (see Table 3). Thus, 
reducing their uncertainties further would not significantly 
affect the calculated D/H uncertainty. Therefore, the main 
source of D/H uncertainty originates from the d(p, 7)? He re- 
action rate. This conclusion is consistent with the findings 
of the recent LUNA studies [29]. Here, the “Original” line 
(i.e., M.F = 1 case) shown indicate the results calculated with 
the default rates and their uncertainties adopted in the original 
PRIMAT code [14], i.e., corresponding to a D/H uncertainty 
value of 1.83%. 


In addition, we found that changing the uncertainty of 
the reaction rate also affected the central value of deuterium 
abundance. We believe that this is because the sampling of 
the reaction rates in the code was performed according to the 
log-normal distribution. This distribution is a typical skewed 
distribution, where the larger the parameter ø, the more the 
most likely position is to be skewed to the left. We found 
that the change in the central value of the deuterium abun- 
dance caused by increasing uncertainty was equivalent to the 
change in the central value caused by the shift in the most 
likely position associated with uncertainty. 


8 mÁ 
—=— p(n, y)d 1.8 
ze dle, He} 5. 
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Fig. 1. Calculated D/H uncertainties v.s. multiplying factor of the 
reaction-rate uncertainty for the four most important reactions. The 
connecting lines are used just to guide the eyes. Table 4 can be ref- 
erenced for details. The inset shows the enlarged plot for the multi- 
plicative factors of M.F. = 0.1, 0.8 and 0.9. Here, the “Original” line 
(corresponding to a D/H uncertainty of 1.83%) indicates the results 
calculated with all the default rates and the associated uncertainties 
adopted in the PRIMAT code [14]. 
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Table 4. Calculated primordial deuterium abundance and uncertainty (in the parenthesis) with different multiplying factors. That is to multiply 
the corresponding default reaction-rate uncertainty by factors of M.F. = 0.1, 0.8, 0.9, 1, 2, 4 and 6, respectively. Here, we put more digits for 
M.F. = 0.1, 0.8 and 0.9 columns for readers to see the tiny differences. 


D/H mean x 107° (uncer. in %) 


M.F=0.1 M.F.=0.8 M.F.=0.9 M.F.=1 (Original) M.F.=2 M.F.=4 M.F.=6 

p(n, y)d 2.4585(1.8234%) 2.4586(1.8250%) 2.4586(1.8255%) 2.459(1.83%) 2.459(1.83%) 2.459(1.86%) 2.459(1.90%) 
d(p,y)?He  2.4590(1.4073%) 2.4587(1.6832%) 2.4586(1.7520%)  2.459(1.83%)  2.457(2.73%) 2.453(4.78%) 2.447(6.83%) 
d(d,n)*He  2.4581(1.7561%) 2.4585(1.7973%) 2.4585(1.8107%)  2.459(1.83%)  2.459(2.06%) 2.460(2.83%) 2.461(3.77%) 
d(d, p)t 2.4584(1.7685%) 2.4585(1.8038%) 2.4586(1.8142%) 2.459(1.83%) 2.459(2.01%) 2.459(2.61%) 2.459(3.36%) 

3 He(n, p)t 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.82%) 2.459(1.83%) 2.459(1.83%) 
t(d,n)’He 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%)  2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 
3He(d,p)'He  2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%)  2.459(1.83%) 2.459(1.84%) 2.459(1.85%) 
3He(a, y) Be  2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 
tla, y) Li 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


2.459(1.83%) 
2.459(1.83%) 


IV. ACCESS TO THE OBSERVATIONAL PRECISION 


In recent astronomical observations [10], the recom- 
mended value of the deuterium abundance is D/H = 
(2.527+0.030)x10—° with an accuracy of approximately 
1.19%, which is better than the above PRIMAT calculated 
accuracy of 1.83%, that is, (2.459+0.045) x 107%. Therefore, 
we focus here on the four reactions that have the largest in- 
fluence on the D/H uncertainty and use a binary search to 
confirm whether reducing their rate uncertainties can improve 
the BBN calculated accuracy to the level of the observation. 
However, we found that the accuracy of the current deuterium 
abundance observations could not be achieved by reducing 
only the uncertainties of the reaction rates, and the uncertainty 
of 0h? should also be reduced. Note that the d(p, y)?He 
reaction rate is adopted from Ref. [22] for all the abovemen- 
tioned calculations. 


Most recently, Mossa et al. [29, 30] measured the 
d(p,y)? He cross section in the E..m, = 33-263 keV energy 
region using the LUNA 400 kV accelerator to an unprece- 
dented precision of better than 3% by exploiting the million- 
fold reduction in cosmic-ray muons at Gran Sasso. Their new 
astrophysical S factor remarkably improved the evaluation 
of the present-day baryon density, 0,h?, using the standard 
BBN model alone. In this study, we utilized the updated re- 
action rate of d(p, y)?He from LUNA [29] to calculate the 
abundances and uncertainties of these primordial light nu- 
clides. The results are listed in Table 5, where the abundance 
of tHe is expressed using its mass fraction Yp, whereas the 
abundances of the other nuclides are expressed as the ratio of 
their number density to that of 'H. We find that when utilizing 
the reaction rate of d(p, y)? He reported by LUNA, the calcu- 
lated D/H uncertainty can be reduced from 1.83% to 1.56%, 
but it does not reach observational precision. We found that 
this goal can be achieved by reducing the current uncertainty 
in Qyh? by approximately 50%. 
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Table 5. Abundances of primordial nuclides calculated with the 
recent LUNA d(p, y)?He rate [29]. The observational values are 
adopted from Ref. [14]. 


Observation (uncer.) Present calc. (uncer.) 


Yp 0.2449(1.633%) 0.2471(0.075%) 
D/Hx1075 2.527(1.19%) 2.471(1.56%) 
3He/Hx10-°  <1.1(18%) 1.0697(2.12%) 


TLVHx1071°  1.58(17.7%-22.2%) 5.5676(4.24%) 


V. CONCLUSION & OUTLOOK 


In this study, we investigated the primordial deuterium 
abundance and its uncertainty using the BBN code PRIMAT. 
We found that the predicted deuterium abundance uncertain- 
ties were dominated by the reaction rate uncertainties in the 
four most important reactions of d(p,y)?He, d(d,n)*He, 
d(d,p)t, and p(n,)d, as well as that in Q,h?. Although 
the current BBN calculation can reach a deuterium abun- 
dance precision level of 1.56% with the recent precise LUNA 
d(p, y)? He rate, there is still a gap of 0.4 % in the observa- 
tional precision. We found that this gap cannot be largely re- 
duced by only reducing the uncertainties of the reaction rates 
and the uncertainty of Qh? should also be reduced. If the 
uncertainty of Q,h? adopted from the Planck 2015 results [4] 
is reduced by approximately 50%, the calculated D/H uncer- 
tainty can reach the observational level. 

In nuclear physics, the reaction rate uncertainties for the 
remaining three reactions must be significantly reduced, ex- 
cept for the d(p, y)?’ He reaction. For instance, relevant re- 
actions can be measured directly at the China Jinping Un- 
derground Laboratory (CJPL) [31-33], which is the deep- 
est operational underground laboratory for particle and nu- 
clear physics experiments worldwide. With such a unique su- 
perlow background environment[34], several successful ex- 
perimental campaigns [35-42] were conducted at the Jin- 
ping Underground Nuclear Astrophysics Experimental Facil- 
ity (JUNA) [43, 44]. 

Figure 2 shows the observational and BBN model calcu- 
lated primordial deuterium abundances and their uncertain- 


257 ties. The observed deuterium abundance was recommended 
as D/H = (2.527+0.030)x10~° [10, 14], whereas the cur- 
rently calculated value was (2.471-40.039)x10~° and there 
was an approximately 1.50 deviation in the central abundance 
value. In addition, the calculated values were less accurate 
than the observed values. Therefore, From a nuclear physics 
perspective, these important reactions still need to be mea- 
sured precisely. Furthermore, the BBN calculated precision 
25 for other primordial nuclides can also be effectively reduced 
2 using more precise reaction rates, with which the cosmic pa- 
2e7 rameters can be strictly constrained. 


258 


259 
260 
261 
262 
263 
264 


3.4 
-—@ Observation 
-—@ Calculation 
3.2 
3.0 


D/H (x10-5) 


Pil 


fF PY gh BO A = & Se & Ss 
< RS . C Ry oe oO N) 
C K) N SP eh FF SI WD > O V2) 
Coe SoH S O O O O O R o ow ew cos 


Fig. 2. Observational and model calculated primordial deuterium 
abundances and the associated uncertainties. The data (from left to 
right) are adopted from Refs. [7, 10-12, 14, 45-53]. The current 
result shown is calculated by using the recent LUNA d(p, yH e 
rate, and by considering the uncertainties in Qh? and Tn. 
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MAT code adopts only a very small (0.45-0.73)% reaction 
rate uncertainty for this p(n,7)d reaction based on theoreti- 
cal calculations [21]. Such a small uncertainty requires fur- 
ther validation using more precise experimental data in the 
BBN energy region. We are now planning to conduct pho- 
todisintegration measurements at the Shanghai Laser Elec- 
tron Gamma Source (SLEGS) facility [63, 64], which uses 
the back-Compton scattering of electrons from the 3.5 GeV 
electron storage ring of the Shanghai Synchrotron Radiation 
Facility (SSRF) and the CO% laser to generate the y beam in 
the energy range of 0.4~20 MeV [65]. The y flux is expected 
to reach 10° ~ 107/s, depending on both the y-ray energy 
and the collimator size [66, 67]. 

In the proposed experiment, a 47 flat-efficiency 3He neu- 
tron detector array was intended for use as a photoneutron de- 
tector with a detection efficiency of approximately 42% [68] 
in the energy region of interest. An LaBrg scintillator detec- 
tor was used for y flux monitoring and waveform acquisition. 
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0.01—10 MeV (as shown in the shaded region in Fig. 3) with 
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abundance uncertainty associated with the p(n, y)d rate. As 
stated in Ref. [21], “the results of the R-matrix theory [69] at 
E =0.1 and 1 MeV differ significantly from the other theoret- 
ical estimations by ~4.6%. Therefore, it would be important 
to experimentally measure the np-capture cross sections at 
these energies to resolve this significant discrepancy”. Thus, 
precise data can also help constrain nuclear reaction models. 

Finally, it is observed that different BBN codes generated 
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different nuclear-physics inputs used. Therefore, the con- 
sistencies of different BBN codes can be verified using the 
same inputs. In addition, the Monte Carlo-type PRIMAT code 
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hence, the Monte Carlo method should be implemented in 
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tant not only to check how the uncertainties in the reaction 
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central rate values may influence the D/H values. All these 
issues, which are beyond the scope of this study, need to be 
studied in the future. 
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